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Foreword 
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00399) and conditionally accepted in July, 1998 pending submission of a minor revision. 
This was withheld by the author because it was simultaneously submitted as part of a 
patent disclosure, LAR-15827, in April, 1998. 

Much has changed since the time this report was written. This model has since been 
improved upon over the past 9 years since it was written. It is being reproduced here for 
historical purposes and to make it more referenceable for those who may wish to 
reference it. 



Blocked Force and Loading Calculations for LaRC Thunder Actuators 


Joel Campbell 

NASA Langley Research Center 
Hampton, VA 23681 


Abstract 

An analytic approach is developed to predict the performance of LaRC Thunder 
actuators under load and under blocked conditions. The problem is treated with the Von 
Karman non-linear analysis combined with a simple Raleigh-Ritz calculation. From this, 
shape and displacement under load combined with voltage are calculated. A method is 
found to calculate the blocked force vs voltage and spring force vs distance. It is found 
that under certain conditions, the blocked force and displacement is almost linear with 
voltage. It is also found that the spring force is multivalued and has at least one 
bifurcation point. This bifurcation point is where the device collapses under load and 
locks to a different bending solution. This occurs at a particular critical load. It is shown 
this other bending solution has a reduced amplitude and is proportional to the original 
amplitude times the square of the aspect ratio. 

Introduction 

Thunder actuators are devices constructed from an isotropic laminate of aluminum, 
LaRC SI adhesive, piezo-electric PZT, and a metal backing of either steel or brass. The 
LaRC adhesive is solid thermoplastic and the aluminum layer is indented to provide 
electrical conductivity after the actuator is cured at 250 deg C. After the curing process 
they form a bow shape after reaching room temperature. If simply supported or arranged 
in a clamshell arrangement they produce a linear motion when an external voltage is 
applied. 

A simple method for predicting the performance of these devices has long been 
sought. As of the time of this writing (2/98), a method for predicting the performance of 
these devices under loading or blocked conditions has never been attempted. One method 
that is widely used to measure blocked force is to place a mass on top and measure the 
displacement under load. This is not true blocked force, however. Blocked force is the 
change in force when constrained not to move. In particular, displacement under load is 
not the same as the change in force when constrained not to move. Because of the 
peculiar non-linear properties of the device, it is easy to be lulled into thinking one is the 
other when, in fact, they are different things. 
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If the device were perfectly linear, the delta amplitude with voltage would be 
completely independent of load in the static limit. This is easy to see by taking a linear 
spring and applying a static external force to find the change in amplitude due to that 
external force under different loading conditions. A non-linear device, on the other hand, 
is quite different and does depend on load. In this case there is actually some critical load 
where the device switches from one bending solution to another at a far reduced 
amplitude. The current blocked force procedure is to apply enough load until the 
amplitude is small. This procedure results in a gross overestimate for the blocked force 
because it depends more on the structural properties of the device than the true actuation 
force. In fact, the true blocking force has never been measured for these devices. 

The correct method to measure blocked force is to constrain the device such that the 
bottom is constrained to move in one plane but fixed to that plane and the top is 
connected to a load cell. In that situation the force will be F=F(x,V). This will be a 
combination of the spring force plus the actuation force. In a particular limit this may be 
approximated by F(x,V)=-k (x-x())+o. V. Of course, this neglects hysteresis effects. 


The incorrect method is to place a mass on top and measure the delta amplitude with 
respect to voltage. This method represents more what it can support structurally. It is a 
common method with these devices to place a mass on top, then apply a low frequency 
sinusoid and measure the delta amplitude. This is usually interpreted as the actuation 
force being the same as the weight of the mass. In reality, this is anything but the case. 

One approach to this problem is to use classical lamination theory [1]. However, this 
leads to an incorrect qualitative behavior [3], As a result, a non-linear approach must be 
used. The type of correction used here is one that is valid for large displacement and 
small strain. This correction is usually attributed to Von Karman [2] in plate theory. Hyer 
[4] has applied this to his treatment of orthotropic laminates with good success. Campbell 
[3] extended this to Thunder actuators by treating each layer as isotropic and including 
voltage terms. 

Solution by the Von Karman Non-Linear Approximation 


The Von Karman approximation assumes large displacement and small strain. This 
suggests including second order terms in w in addition to the linear terms for u, v and w. 
Expressing the strain to include these terms gives, 


~ llli , J_ 
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By using the relation, 
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The total energy takes the form, 

U= j jj ( dU 0 -dW T -dW II w p(x, y) dx dy, (4) 


where dU l is the stored elastic energy volume density and dWj and dWA are the 
temperature and actuation contributions. The last integral is the work done against the 
external pressure, p. The elastic energy portion of it for the k'th layer is, 



(5) 


For the thermal portion we have, 

(dW% k = (a T x e x +o T y e) 


(<2i + Q^aAT s v + (<2 1 + <2 2 ) a AT e v 


( 6 ) 


For the actuation portion we have, 


3 



( 7 ) 


= (°v e x + o y e^ _ 

= (<2 1 + Q 2 ) A\ Y £ x + (01 + 2 2 ) ^31 y S v 

Actuation is assumed to occur on the j'th layer and Q is the reduced stiffness with 
Ql l=Q22=Ql and Ql2=Q21=Q2- Now that we have the general form for the energy 
integral we may try a Raleigh-Ritz type solution. The method involves making a guess 
for the solution. This guess should look qualitatively like the expected solution yet be 
flexible enough to allow for adjustment. This adjustment comes in the form of changing 
the guess based on a finite number of parameters which the guess contains. The 
parameters are adjusted in such a way as to minimize the total energy of the system. In 
this situation it would be advantageous to choose a guess that is close enough to the 
classical lamination result that in the linear limit (small scaling) the solution approaches 
the classical lamination result. Such a guess takes the form, 

w° =w% + ^(a x 2 + b (8) 

Classical lamination theory assumes no shearing strain between layers for thermal 
expansion. If we make the same assumptions here we find u^ and v^ must take the form, 

u° = c x- ^ a 2 y? - ^ ab r/ 

(9) 

v° = d y - -i b 2 f - i a b x 2 y 

The ansatz represented by Equation 27 and 28 has been used successfully by Hyer [3] 
in his treatment of orthotropic laminates and by Campbell [3] with Thunder actuators. It 
is sufficiently flexible enough to be used with isotropic laminates as well as we will soon 
see. Equations 3, 8 and 9 give, 

e x =c~^abf-az 

z y =d-^abx 2 -bz (10) 


If the length in the x direction is L x and the length in the y is Ly then, we may 
approximate w 0 ° by. 
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(11) 


<=-g-(« ^r + b ity 

This assumes the device is simply supported. For p(x,y) we choose a point load of, 

p(x, y) =F 8 (x)6 (y) (12) 

Other choices are possible such as a constant distributed pressure. However, a point 
load is closest to the type of load used with the rectangular devices. The solution is 
obtained by minimizing the total energy with respect to a, b, c and d. The result of the 
minimization is, 

D x a + K 8 b - B x c - B 1 d + K x b 2 + K 6 ab - K 4 b c - K^b d + 

K, a b 2 + M' - i F= 0 

K & a + D x b - B 2 c - B x d + K 2 a 2 + K 5 ab-K 4 ac - K 2 a d + 

F = 0 (13) 

- B x a - B 2 b + A 1 c+A 2 d-K 4 ab-N' = 0 

- B 2 a - B x b + A 2 c + A x d - K 2 a b - FT = 0, 

where, 
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The A's, B's and D's have the usual meanings of extensional, coupling and bending 
stiffness. The force and moment per unit width, N' and M', are given by, 
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where, 
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and. 


N T = AT Jj a, (<2, + Cz), [z t - 4. ,) 

^ e.+sD. Gs-^-D 


(17) 


The z's are measured with respect to the middle surface and zo is the top surface. The 
lengths L x and Ly are the length and width of the actuator. It is a simple matter to show 
that the solution in the linear limit (small scaling - L x =Ly=0) matches the classical 
lamination result if we take a=-K x , b=-Ky, c=e x and d=Sy and F=0. 

Unlike the classical lamination result, the non-linear analysis predicts unequal 
curvatures that depend on the magnitude of the scaling and the aspect ratio. It also 
predicts multiple possible solutions. The nature of these solutions are dome-like solutions 
which approach cylinders in the limit of large scaling. Under load, saddles are also 
possible. Typically, there will be one or two stable solutions and possibly a third unstable 
solution depending on the size and material properties. There are also unphysical 
solutions in terms of complex numbers. 

As a simple example we take 3 layer actuator constructed from 1 mil brass, 1 mil 
LaRC Si thermoplastic and 6.8 mil PZT 5A. By solving Equation 13 with V=0 and 
plotting F as a function of dome height we find the spring force for the actuator. The 
results of this calculation show the spring force is multivalued. When initially in the short 
axis mode, the actuator may be depressed until a critical load is reached. At this point, the 
actuator switches and locks to the long axis solution. This has been verified 
experimentally and is a general characteristic of all Thunder actuators that have unequal x 
and y curvatures. When the x and y curvatures are strongly unequal, the spring force is 
almost linear except for the bifurcation point. 
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Figure 1 . Calculation of spring force for a 3 layer actuator. 


Blocked force as a function of voltage may be calculated by constraining the distance 
by Equation 11. and solving Equation 13 to find the resulting force as a function of 
voltage at a particular constrained distance. Figure 2. is the result of such a calculation for 
the same 3 layer actuator. In this case the blocked force appears to be almost linear with 
distance. In reality there is some hysteresis involved. The horizontal axis is voltage and 
the vertical is pounds. The constrained distance chosen was close to the equilibrium point 
and the short axis solution is shown in the plot. 



Blocked Force vs. Voltage 

Figure 2. Calculation of blocked force as a function of voltage for a 3 layer actuator. 

If the constrained distance is close enough to the bifurcation point it is actually 
possible to cause the actuator to switch from the short axis solution to the long axis 
solution. The result is a bifurcation point involving voltage and blocked force. 

Special Limits 

There are two special limits that reduce to analytic forms. The first is the classical 
lamination limit where the size to thickness ratio is very small. The second is the strongly 
non-linear limit where the x and y curvatures are strongly unequal. This second limit is 
the most common type of rectangular actuator encountered. However, other types are 
possible and it is impossible to say for sure how the actuator will turn out without first 
doing the non-linear calculation first. It is highly desirable to design an actuator that is as 
close to a cylinder as possible. This is because it will behave more predictably than one 
that is close to the curvature vs. size bifurcation point. Actuators that are designed to be 
valid in the classical lamination theory limit are predictable but have a much reduced 
amplitude. 
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To find the classical lamination theory limit we solve Equation 13 in the limit Lx=0, 
Lx=0 with Lx/Ly fixed. The result is, 


(a + a) | 

AT -1 

^)fJ-(F i+j B) W 

(A + A) 

(a + a) 

-(s,+sy 


(A + At) | 

f , 

Af- 

fJ - (fl, + fi) /V 

(A + A 2 ) 

(A + A) 

- ( B , + B$ 


(18) 


Using the results of Equation 18 with the results of Equations 15, 16 and 17 to find F 
expressed in terms of the dome height, as, 


where, 


F=-£(!-i 0 ) + p v 


(19) 
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The other limit is the strongly non-linear, perfect cylinder limit. In this case Equation 13 
becomes independent of either a or b, giving two possibilities. These two possibilities are, 


D x a - B x c - B 2 d + Af - [ k] F= 0, 


— Fj r/ + A| c + At d — N' — 0, 


( 21 ) 


— B~, g + At c + A| d — — 0, 


or, 
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(L\ 

D x b - B 2 c - B x d + M \ F ={) 


B-, b + A x c + At d- BF = 0 


( 22 ) 


- B x b + A 2 c + A x d-N' = 0 , 

Equations 21 and 22 may be solved for a or b so that we may find the force as before. 
The result is, 


where, 
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With these results we may estimate the critical load represented in Figure 1. To find 
this we use Equation 23 with V=0 and write Fa=Fb so that, 


k a (§ “ §«) = k b ($ - 3 /) 


(26) 


or, 


k„ - ku 


(27) 


By substituting this into Equation 23 we find the critical load is then, 

F = F = - — — fH - H 1 

a b k -k, 


(28) 


We may also find the static amplitude with voltage under load. There are two 
possibilities, 


mg=-k a (|-^) + fi a V, 
mg = -k b + fi b V 


The amplitudes are, 


- _ £ mg $ 

51 k k ’ 


a a 


= r>-?L* + hv 


(29) 


(30) 


If the device is originally in the first bending solution and the critical load is applied, 
it will switch to the second bending solution. The ratio of the delta amplitude with 
voltage is. 



— , r . 2 
kb fib k a _ f ^ 

fijL A V $ a k b 


(31) 


This reduction in static amplitude with critical load is a purely a structural effect. The 
true actuation force, on the other hand, is given by Equation 23 - not Equation 28. They 
are completely different things. 
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